Stochastic finite differences for elliptic diffusion equations in stratified domains

We describe Monte Carlo algorithms to solve elliptic partial differential equations with piecewise constant diffusion coefficients and general boundary conditions including Robin and transmission conditions as well as a damping term. The treatment of the boundary conditions is done via stochastic finite differences techniques which possess an higher order than the usual methods. The simulation of Brownian paths inside the domain relies on variations around the walk on spheres method with or without killing. We check numerically the efficiency of our algorithms on various examples of diffusion equations illustrating each of the new techniques introduced here.